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We investigate an intra-cavity coupled down-conversion scheme to generate quadripartite entan- 
glement using concurrently resonant nonlinearities. We verify that quadripartite entanglement is 
present in this system by calculating the output fluctuation spectra and then considering viola- 
tions of optimized inequalities of the van Loock-Furusawa type. The entanglement characteristics 
both above and below the oscillation threshold are considered. We also present analytic solutions 
for the quadrature operators and the van Loock-Furusawa correlations in the undepleted pump 
approximation. 
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I. INTRODUCTION 

Entanglement is a concept of central importance in 
quantum theory and continues to inspire both theoreti- 
cal and experimental efforts to explore quantum systems. 
In addition to this, entanglement is the main resource 
of quantum information and in particular, multipartite 
continuous- variable (CV) entangled states have grown to 
be pivotal in multipartite quantum communication ■ 
The criteria which must be satisfied to establish whether 
bipartite entanglement exists in a given system are well 
known for the CV case @, Furthermore, bipartite 
entanglement can be realised experimentally. The crite- 
ria for the bipartite scenario have been generalized for 
multipartite entanglement by van Loock and Furusawa 
fToj l . Advances have also taken place in the experimen- 
tal generation of tripartite entanglement. In particular, 
there have been experiments where entangled beams are 
produced by mixing squeezed beams with linear optics 
IS, HL IE S EH El ■ Several proposals have also been made 
whereby multi-frequency entangled outputs are gener- 
ated. These rely on the use of non-degenerate down- 
conversion |13| or cascaded or concurrent nonlinear opti- 
cal processes |14| - fl8j where the tripartite entanglement is 
instead produced via the interaction with the nonlinear 
medium itself. It is these concurrent processes that we 
consider in this work for the quadripartite case. 

In regards to quadripartite entangled beams there have 

been theoretical proposals based on linear optics and cas- 
i — 1 1 — ~i 

cadcd nonlincarites [19H22j . In this work we build on a 
tripartite scheme proposed by Bradley et al. in [l6[ but 
for the case of quadruply concurrent nonlinearities. Fur- 
thermore, we use an optimized version of the van Loock- 
Furusawa inequalities to demonstrate quadripartite en- 
tanglement in this system. 

This paper is organized as follows. In Section II we 
describe the Hamiltonian for the system under consid- 
eration. Section III discusses the van Loock-Furusawa 



(VLF) criteria as a means of quantifying quadripar- 
tite entanglement. Section IV considers the interaction 
Hamiltonian in the undepleted pump approximation and 
gives analytic solutions for the quadrature operators, as 
well as the VLF correlations. In Section V we present the 
full equations of motion for the system and Section VI 
gives the steady state solutions to the classical versions of 
the equations of motion and provides an overview of the 
linearized fluctuation analysis used in this work to calcu- 
late the measurable output fluctuation spectra from the 
cavity These output spectra are are found in Section VII 
and used to demonstrate violation of the optimized van 
Loock-Furusawa criteria and hence, demonstrate quadri- 
partite entanglement. 



II. SYSTEM AND HAMILTONIAN 

We model a system in which pump lasers drive four 
modes in an optical cavity. As depicted in the simpli- 
fied experimental setup shown in Figure [TJ the four in- 
puts interact with a \^ non-linear crystal to produce 
four low-frequency entangled output beams at frequen- 
cies u>5, uje, W7, u>s- For example, mode 1 is pumped at 
a particular frequency and polarization such that it pro- 
duces modes 5 and 6. 

The full Hamiltonian for the eight-mode system, de- 
scribing the interaction inside the optical cavity and the 
interaction of the cavity fields with the output fields, can 
be written as 

^ 'rtpump ~t~ 7~Lint "F 7~Lfree "F 7~Lres-> (1) 

where the interaction Hamiltonian is 

Hint = ih[xiaia\a\ + X2a2'a\a\ + x^^ol 
44] 



XAaA<i\a\] + H.c, 



(2) 



with the Xi representing the effective nonlinearities and 
hi denoting the bosonic annihilation operators for the in- 
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FIG. 1: (Color online) Schematic of a x crystal inside a 
pumped Fabry-Perot cavity. Pump lasers drive four intracav- 
ity modes with frequencies wi, 0J2, ^3 and W4 (represented by 
circles and squares) , which are down-converted to four output 
modes with frequencies cj 5 , ujh, 0J7 and cjg. 



tracavity modes at frequencies u>i. The pumping Hamil- 
tonian, describing the the cavity driving fields, in the 
appropriate rotating frame is 



Hpump ^ 



e { cn 



i=l 



(3) 



and the cavity damping Hamiltonian is given by 



8 



[fid] + f l&i 



(4) 



where are the classical pumping laser amplitudes for 
mode i, and the Ti are the annihilation operators for 
bath quanta to which each of the intracavity modes are 
coupled, and which represent losses through the cavity 
mirror. 



Using the quadrature definitions in Equation ([5]), the 
quadripartite inequalities which must be simultaneously 
violated by the low frequency modes are, 

V(X 5 - X 6 ) + V{% +%+ g 7 Y 7 + g s %) > 4, (6) 
V(X 6 - X 7 ) + V{g b % +%+%+ 98%) > 4, (7) 
V(X 7 - X 8 ) + V{g 5 % + g 6 % + Y 7 + %) > 4, (8) 

where the gi(i = 5,6,7,8) are arbitrary real parameters 
that are used to optimize the violation of these inequal- 
ities. In particular, we minimize Equations © and (|5]) 
with respect to g 7t8 and g$fi, respectively. Solving the 
resulting equations leads to the optimized expressions, 
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(9) 
(10) 

(11) 
(12) 



where for covariances we use the notation Vij = ((YiYj) + 
(YjYi})/2 — (Yi)(Yj) and for the case where i = j the 
covariance, denoted Vi, reduces to the usual variance, 
V(Yi). It is important to note that in the uncorrelatcd 
limit these optimized van Loock-Furusawa criteria ap- 
proach 4. Hence, without optimization some entangle- 
ment which is present may be missed. 



IV. ANALYTIC SOLUTIONS IN THE 
UNDEPLETED PUMP APPROXIMATION 



III. QUADRIPARTITE ENTANGLEMENT 
MEASURES 



In order to investigate multipartite entanglement, and 
in particular show that the system under consideration 
demonstrates true quadripartite entanglement, we first 
define quadrature operators [3(| for each mode as, 



Xi 



Y, 



-i(a,i - aj) 



(5) 



such that [Xi, Yj] = 2i. Based on this definition V{Xi) < 
1, for example, indicates single- mode squeezing where 
V{A) = (A 2 ) — (A) 2 denotes the variance. The condi- 
tions proposed by van Loock and Furusawa [Toj . which 
are a generalization of the conditions for bipartite entan- 
glement [1, Q are sufficient to demonstrate multipartite 
entanglement. We now demonstrate how these may be 
optimized for the verification of genuine quadripartite en- 
tanglement in this system. 



It is useful to consider the interaction Hamiltonian in 
the undcplctcd pump approximation in the absence of 
a cavity in advance of a more complete approach that 
considers the full quantum equations of motion for all of 
the interacting fields inside a cavity. We stress here that 
these equations are not of exact physical relevance, but 
do give useful insights into the properties of the Hamilto- 
nian. Here we show that it is possible to obtain analytic 
solutions for the quadrature operator equations of motion 
using the undcplctcd pump approximation. This entails 
setting £i = Xi(ai(0)) (i=l,2,3,4), where & are positive, 
real constants. Under these conditions the interaction 
Hamiltonian can be written as, 



Hi 



ih 



f At 



+ &(a 7 % - a 7 a 8 ) + &(a T 8 al - a 8 a 5 ) . (13) 



The Heisenbcrg equations of motion can then be writ- 
ten, 
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where 



dag 
~~dt 

~dt 

dd 7 

~dt 
das 
~dt 



and these equations can be recast in terms of the quadra- 
ture operators as follows, 



dX* 



= 64 + 64' 


(14) 






= 64 + 64> 


(15) 


.4 


= cosh 2 (^), 


= 64 


(16) 


B 


= isinh(2^) 


= 64 + 64> 


(17) 


C 


= sinh 2 (£i). 



dt 
d% 

dt 
dX 6 

dt 
d% 

dt 
dX 7 

dt 
dY 7 
~~dt 
dX 8 

dt 

dYs 
dt 



6*6 + 


■ 6*8, 


(18) 


-6*6 


" 6^8, 


(19) 


6*5 4 


" &X7, 


(20) 


-6*5 


~ 6*7, 


(21) 


6*6 4 


" £3X3, 


(22) 


~6*6 


- 6^8, 


(23) 


6*7 4 


- 6*5, 


(24) 


-6*7 


- 6*5- 


(25) 



It is these equations that we solve to find analytic solu- 
tions for the quadrature operators as functions of their 
initial values. 



A. Solutions with equal £j 

To begin with, we set all the interactions equal so that 
£i = £ and find analytic expressions for the VLF correla- 
tions by solving the Heiscnberg equations of motion for 
this case. The solutions for the quadrature operators are 
found to be, 



(34) 
(35) 
(36) 



From these expressions for the quadrature operators it 
is possible to find the variances and covariances neces- 
sary to calculate the van Loock-Furusawa (VLF) criteria 
within the undcplctcd pump approximation. In fact, the 
variances are all equal and given by the following time- 
dependent moments, 



(Xf) = (Y 2 ) = A 2 + 2B 2 + C 2 , 



(37) 



since the expectation values of the amplitudes are all 
zero. Here we have used the fact that (Xi(0)Xj(0)) = 
(Yi (0)Yj (0)) = Sij . A similar approach can be used to cal- 
culate the covariances which are equivalent to the time- 
dependent moments (XiXj) and (YiYj). In particular, 
the covariances are given by, 



(X 5 X 6 ) = -(%Y 6 ) = 2{AB + BC), (38) 

(X 5 X 7 ) = (Y 5 Y 7 ) = 2(AC + B 2 ), (39) 

(A 5 A 8 > = -(Y 5 Y 8 ) = 2(AB + BC), (40) 

(X 6 X 7 ) = -(Y 6 Y 7 ) =2(AB + BC), (41) 

(A 6 A 8 ) = (Y 6 Y 8 ) =2(AC + B 2 ), (42) 

(A 7 A 8 ) = -(Y 7 Y 8 ) =2(AB + BC). (43) 



x 5 (t) 


= AX 5 {0)-i 


- BX 6 (0) - 


t- CX 7 (0) 


4-5X 8 (0)(26) 


Mt) 


= AY 5 (0)- 


B%(0) + 


CY 7 (0)- 


B% (0), (27) 


x 6 (t) 


= 5X5(0)4 


-AX 6 (0)- 


lBXr(0) 


4-CA 8 (0){28) 


Mt) 


= -5*5(0) 


+ a%(o) 


- bMo) 


4-CT 8 (0),(29) 


X 7 (t) 


= CA 5 (0)4 


-BX 6 (0)- 


f AX 7 (0) 


4-5A 8 (0){30) 


Mt) 


= cy 5 (o)- 


B%(0) + 


AY 7 {Q) - 


5*8(0), (31) 


Xs(t) 


= B X 5 (0)4 


-CX 6 {0)- 


f 5A 7 (0) 


4-AA 8 (0)(32) 


Mt) 


= -BY 5 (0) 


+ CY 6 (0) 


- B*V(0) 


+ Ay 8 (0),(33) 



From these variances and covariances we obtain Equa- 
tion @3J an analytic expression for the optimized VLF 
correlations defined in Equations [6][Sj All three VLF cor- 
relations are equal when £j = £ and hence we label any 
of the correlations in Equations [fj][5] as V3. Figure [2] pro- 
vides a plot of these optimized VLF correlations, V3, as 
a function of £t. We observe that quadripartite entan- 
glement is present in this system with V3 < 4 for all £t. 
This suggests a more complete treatment incorporating 
depletion of the pump fields and a cavity, will find where 
quadripartite entanglement is present. 
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V 3 = AA 2 - AAV2B + 4(V2B - C)C 



2{B 2 - 4B 3 / 2 CV2 + 12BC 2 - 8C 3 V2B + AC 4 
A 2 - 2Ay/B + B- G'VZB + C 2 



(44) 




FIG. 2: Analytic solutions for the optimized van Loock- 
Furusawa correlations, V3, found by solving the Heisenberg 
equations of motion in the undepleted pump approximation. 
A value of less than 4 signifies quadripartite entanglement. 
All quantities depicted here, and in subsequent graphs, are 
dimensionless. 



D 

E 

F 
G 
H 
I 
J 



cosh 2 (f2i/%/2), 
sinh 2 (fit/V2), 

£2 sinh(V2m) 

Q + £ 2 cosh(V2^t) 

a6(cosh(72m) - 1) 
ft 2 

tj + £f co&h{y/2Q.t) 



(53) 
(54) 

(55) 
(56) 
(57) 
(58) 
(59) 



We note here the generality of the solutions presented and 
that other cases are possible numerically. The variances 
are given by, 



(60) 
(61) 
(62) 
(63) 



(xl) 
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(xl) 
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B. Solutions with unequal £4 



and the covariances are, 



For simplicity, here we assume that £1 = £2 and £3 = £4 
and setting = ysi+Cf we nn d that the solutions are 
given by, 



x 5 (t) 

X 6 (t) 

%(t) 
X?(t) 

x s (t) 



DX 5 (0) + EX e (0) + FX 7 (0) + GV 8 (0)(45) 
D% (0) - E% (0) + FY 7 (0) - G% (0) , (46) 
EX 5 (0) + HX e (0) + EX 7 (0) + 7A 8 (0),(47) 
-EY 5 (0) + H%(0) - EY 7 (0) + I%(fl), (48) 
FX 5 (0) + EX 6 (0) + DX 7 (0) + GA 8 (0)(49) 
FY 5 (0) - E%(0) + DY 7 (0) - GV 8 (0), (50) 
GX 5 (0) + IX 6 (Q) + GX 7 (0) + JA 8 (0), (51) 

-gyb(o) + /y 6 (o) - gy 7 (o) + jr 8 (o), (52) 



(X 5 X 6 ) 
(X 5 X 7 ) 
(X 5 X S ) 
(X 6 X 7 ) 
(X 6 X S ) 
(X 7 X 8 ) 



-(%Y 6 )=DE 
(Y 5 Y 7 ) = 2DF ^ 
-(Y 5 Y S )=DG 
-(%Y 7 ) = DE 
(Y 6 Y 8 ) = 2EG ^ 
-(Y 7 Y S ) =DG 



-EH + EF 
E 2 + G 2 , 
-EI + FG- 
- EH + EF 
HI + IJ, 
-EI + FG- 



f IG, (64) 
(65) 
GJ, (66) 
f IG, (67) 
(68) 
GJ. (69) 



where 



As described in Section llV A[ we can now calculate the 
VLF correlations. Figure[3]shows the optimized VLF cor- 
relations as a function of £ii, for the case £2 = 0.5£i. 
Again, we observe that quadripartite entanglement is 
present in this system with all three VLF correlations 
less than 4 for some £ii. Comparing Figures [5] and [3] we 
see that the greatest degree of entanglement is obtained 
for the case where all £, are equal. Figure [3] shows that 
for each VLF correlation the entanglement is degraded 
beyond a particular value of £ii, however this is not the 
case in Figure [3] when all the £i are equal. 
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FIG. 3: Analytic solutions for the optimized van Loock- 
Furusawa correlations, V56, Vq7 and Wg with £2 = 0.5£i, found 
by solving the Heisenberg equations of motion in the unde- 
pleted pump approximation. Having all three of the correla- 
tions drop below 4 is sufficient to demonstrate quadripartite 
entanglement. 



EQUATIONS OF MOTION FOR THE FULL 
HAMILTONIAN 



where ji are the cavity loss rates at the respective fre- 
quencies and T>i[p] = 2hipa\ — djdj/3 — pd\di is the Lind- 
blad superoperator [23| under the usual zero-temperature 
Markov approximation. From this one can derive the 
stochastic differential equations (SDEs) in the positive- 
P representation [24|, and in turn, study the intracavity 
dynamics. 



Our approach [25| involves converting the quantum op- 
erator equations of motion of Equation [70] into a Fokkcr- 
Planck equation for the positive- P representation pseu- 
doprobability distribution of the system [24], [2f|. This 
can then be interpreted as a set of c- number SDEs. It 
should be noted that the use of the positive- P repre- 
sentation, rather than the Glauber-Sudarshan P repre- 
sentation (2(| is necessary to ensure that the dif- 
fusion matrix of the FPE is positive-definite. This is 
achieved with the positive-P approach by defining two 
independent stochastic fields cti and af corresponding to 
the mode operators dj and dj, respectively, in the limit 
of a large number of stochastic trajectories. Using this 
method it is possible to calculate stochastic trajectory av- 
erages which correspond to the normally-ordered expec- 
tation values of quantum-mechanical operators, for ex- 
ample, (ai) n (a^) m = ((dj) m d™). Taking this approach 
yields a diffusion matrix of the form, 



We will now consider the full physical system, where 
the nonlinear media are contained inside a pumped res- 
onant Fabry-Perot cavity. The master equation for the 
density operator of the system can be found in the stan- 
dard manner by tracing over the reservoirs [23[ and is 
given by, 



D 




d 



(71) 



dp 
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H, 



pump 
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H int ,p\ +Y J 1iT>i[p] (70) 
=1 



where is an 8 x 
is given by, 



null matrix and the non zero block 
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(72) 



The matrix d can be factorised such that the Ito stochas- 
tic differential equations are obtained. For the high fre- 
quency fields this process yields, 



do>i 
~dt 
da 2 
~dT 
das 
dt 
da^ 
~dT 



(2 - X2CtQO> 7 - 72«2, 



£3 - X3«7a8 - 73<*3, 



e 4 - XiCt & a 5 - 74CX4, 



(73) 
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and also the equations found by interchanging a, and a~l . 
Whilst for the low frequency fields one obtains, 



da 5 
~dT 



doe 
dt 



da*? 
~dT 



da s 
dt 



75«5 



(^13(0 + iVu{t)), 



Xi^ia^ 



76"6 



X2«2 



(r] 9 (t) +ir] 10 (t)), 



X3a 3 aj 



77«7 + 
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- «7io(*)), 



X3a 3 a^ 



X4«4a^ 



78"8 



X3"3 



X3a3 



(%(*)+ "76 (*)) 



(%(*) - 



(Vi(t) + i m (t)) 



da 1 
~dt 
da 2 
~dT 
da 3 
dt 
dot4 

~dT 

dct 5 
dt 
da,Q 

~dt 
daj 
~dT 
da 8 
dt 



e l - Xl a 5&6 - 71 0-1, 




£2 - X2<56«7 - 72«2, 




£3 - X3Q7Ci8 ~ 73«3, 




e 4 - X4 a 3&5 ~ 74<5 4 




Vi n'^ -1- v^rv^rv^ — 


15 Lt 5 ? 


Xl"l«5 + X2Q!2<57 - 


76 a 6 , 


X2a 2 a* 6 + X3a3<5g - 


770:7, 


X3a 3 ai + XiOtiol - 


78^8, 



(75) 



+ 



X4"4 



(Vi3(t) - ir)i4{t)), 



and also the equations found by interchanging Oj and 
a+ and ?7,(i) and 77^+2 (t). The 7j are the cavity loss 
rates at the respective frequencies, and r)i{t) are real, 
independent, Gaussian noise terms which satisfy 77,(4) = 
and T]i(t)rjj(t') = SijS(t — t'). It should be mentioned 
that we are assuming that all the intracavity modes are 
resonant with the cavity, and although it is possible to 
include detuning we do not do so here. 



and from these we can obtain steady-state solutions. 
In the remainder of this paper we consider a symmetric 
( r ?i(0 — iV2{t)) system where all the high frequency modes have the same 
cavity damping rates, with 7; = 7 for i = 1,2,3,4 and 
all low frequency modes also have equal cavity damping 
rates, with 7.; = k for i = 5,6, 7, 8. In addition to this, 
(74) we assume all the nonlinearities and hence all the pump 
field amplitudes are equal, that is, Xi = X an d = e, 
respectively. 

For this completely symmetric system, we verify that 
there is an oscillation threshold at the critical pumping 
amplitude, 



"fK 

2X 



(76) 



as is the case for triply concurrent downconversion |16j |. 
This result differs from the standard non-degenerate 
OPO threshold condition by a factor of a half. This 
difference here is due to the fact that each pump mode 
drives two down-conversion processes. The stationary 
solutions below this threshold value are found to be, 



VI. STABILITY ANALYSIS AND 
FLUCTUATION SPECTRA 



e 

a s = - 

7 

&i = 



for % G {1,2,3,4}, 
fori G {5,6,7,8}. 



(77) 



Whilst, above threshold the stationary solutions are 
given by, 



We conduct a linearized fluctuation analysis [23j of the 
system for the purposes of calculating the output spectral 
correlations for the cavity from the intracavity spectra. 
We begin by neglecting the noise terms in Equation (|74|) 
so that af — > a*, and also set a, = + Sa i7 where ocj, is 
a mean value and 5at represents the fluctuations. This 
gives a set of classical equations for the mean values, 



on = — fori G {1,2,3,4} (78) 



Oi = y/(e-e c )/x for »G {5,6,7,8}. (79) 

We see that the low frequency modes become macro- 
scopically occupied as the pumping is increased and the 
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high frequency modes remain at their threshold value. 
Using Equations [73] and [73] we also perform dynamical 
simulations to confirm the steady state values. 

We then proceed to study fluctuations around the 
steady state which allows one to calculate measurable 
output fluctuation spectra plj and hence, quantify the 
quantum correlations of the system. The linearized equa- 
tions for the fluctuations are of the form, 



dSa = -ASadt + BdW, (80) 
where Sa = [5ai, Sd{, 8a%, Sa^ , . . . , Sas, <5aj|~] T , B is the 



noise matrix of Equation [74] with the steady-state values 
inserted, dW is a vector of independent, real Wiener 
increments [25| and A is the drift matrix with the steady- 
state values inserted and given by, 



Ai A 2 

-(a;) t a 3 



(81) 



where Ai = —jhi 
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(83) 



Provided the eigenvalues of the drift matrix A have 
no negative real part the system is stable and we can 
treat the fluctuation equations as describing an Ornstein- 
Uhlenbeck process (28|. This allows one to calculate the 
intracavity spectral correlation matrix, 



S(uj) = (A + iiul)BB T {A T - iujl)- 1 , (84) 



and this is related to the measurable output fluctuation 
spectra using the standard input-output relations for op- 
tical cavities [29j]. Furthermore, it supplies us with all 
that is necessary to calculate the measurable extracavity 
quadripartite entanglement. 



VII. OUTPUT FLUCTUATION SPECTRA 

The same inequalities as given in Section IIIII in terms 
of variances also hold when expressed in terms of the 
output spectra, and it is these quantities that can be 
measured in experiments. In the following, we use the 
notation Zf/*(a;) (ie. any of I§f , Iffl* , Iff) to represent 
the three output spectral correlations of the same form 
as the optimized expressions in Eqns. - ©• In Fig- 
ure 01 we plot these three correlations as a function of 
frequency for our completely symmetric system for the 
below threshold case (solid line) and the above threshold 
case (dashed line). In both cases, the three correlations 
are equal, ie. I§f = I^f = J$". 

The inequalities are violated both below and above the 
threshold, demonstrating quadripartite entanglement. 
Specifically, the results shown are for the case of the 
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pump field amplitudes set at e = 0.8e c and e = 1.2e c . 
For these parameter choices, the largest violation of the 
van Loock-Furusawa entanglement criteria, and thus the 
maximum quadripartite entanglement, is for the low fre- 
quency modes below threshold. In general, for both 
cases the largest degree of violation of the inequalities 
is observed near zero frequency. For large frequencies 
/^"'(w) — > 4 which is the uncorrelated limit for our opti- 
mized expressions. 



5r 

4.5 




co (units of k) 



FIG. 4: The output spectral correlations, /^"'(w), as a func- 
tion of frequency u) (units of n) corresponding to the quadri- 
partite entanglement criteria given in Equation 21 - 23. The 
below threshold (solid line) and above threshold (dashed line) 
cases are shown, for cavity pump amplitudes e = 0.8e c and 
e = 1.2e c , respectively. The three correlations are equal for 
each case for the chosen parameters, which are symmetric 
with \i = X, e i = e ) 7 = 10, k = 1 and \ = 10~ 2 . All 
quantities plotted here and in the subsequent graph are di- 
mensionless. 

We also determine the maximum quadripartite entan- 
glement for the same parameters as in Figure 21 but for 
a range of pump field amplitudes on both sides of the 
oscillation threshold. This is shown in Figure [SJ where 
we plot the minimum value of the output spectra at any 
frequency, as a function of e/e c . As expected [3l|, l32f . 
we observe the maximum quadripartite entanglement at 
the critical pumping amplitude, with the caveat that the 
linearized fluctuation analysis is not valid in the imme- 
diate vicinity of the threshold. The fact that a gradual 
slope is observed below threshold in the region of maxi- 
mum entanglement could prove useful for future experi- 
mental realizations of this scheme. It is also found that 
quadripartite entanglement persists well above threshold, 
with a large violation of the van Loock-Furusawa criteria 
still present as the pumping is increased above the crit- 
ical pumping value. In particular, well above threshold 
the minimum of I ^j ui) approaches 3. This behaviour 
is also seen in |3lL l3^| where the equivalent correlations 
also asymptote to a finite value. 




0.5 1 1.5 2 

e/e 



FIG. 5: Maximum quadripartite entanglement as a function 
of the ratio of the cavity pumping to the pumping thresh- 
old. The cavity parameters are the same as in Figure [4] with 
7 = 10, k = 1 and \ = 10 -2 . Again, all three correlations 
are equal for these parameters. It should be stressed that at 
e/e c = 1 the validity of the results is limited as the linearized 
analysis is no longer valid. 



VIII. CONCLUSIONS AND OUTLOOK 

We have demonstrated intracavity continous-variable 
quadripartite entanglement in quadruply concurrent 
downconversion, both above and below the critical pump- 
ing threshold using optimized van Loock-Furusawa cri- 
teria. Above threshold the proposed scheme produces 
a source of bright entangled output beams. The below 
threshold regime provides the greatest degree of entan- 
glement and a region where this entanglement could be 
measured in experiments. One of the advantages of this 
type of scheme lies in the number of different regimes that 
can be explored by tuning various parameters in exper- 
iments. For example, the pump intensities and coupling 
strengths can be tuned and this makes it possible to vary 
the degree of entanglement in the system. 

Throughout this article, we have studied the properties 
of the interaction Hamiltonian, presented the full quan- 
tum equations of motion and performed a linearized fluc- 
tuation analysis. All results indicate that this system is 
a good candidate for the demonstration of quadripartite 
continuous- variable entanglement. In relation to exper- 
imental implementation of the scheme presented here, 
stabilizing a single cavity in which four entangled modes 
are created may prove preferable in some applications 
rather than alternative schemes which rely on stabilizing 
and synchronizing multiple OPOs. 

Finally, this result could be of further significance in 
a similar system where one of the nonlinear couplings is 
absent. Such a system may be a candidate for realizing 
the simplest four node cluster state [HI, [34| . 
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